Anomalous diffusion of phospholipids and cholesterols in a lipid bilayer and its origins 
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Combining extensive molecular dynamics simulations of lipid bilayer systems of varying chemical 
composition with single-trajectory analyses we systematically elucidate the stochastic nature of the 
lipid motion. We observe subdiffusion over more than four orders of magnitude in time, clearly 
stretching into the sub-microsecond domain. The lipid motion delicately depends on the lipid 
chemistry, the lipid phase, and especially on the presence of cholesterol. We demonstrate that 
fractional Langevin equation motion universally describes the lipid motion in all phases including 
the gel phase, and in the presence of cholesterol. The results underline the relevance of anomalous 
diffusion in lipid bilayers and the strong effects of the membrane composition. 
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Recent advances in single molecule spectroscopy unveil 
anomalous diffusion of microscopic tracers in the crowded 
environment of living cells, starting to reshape our vievifs 
of molecular cell biology and underlining the role of mod- 
ern statistical physics Extensive experimental stud- 
ies show subdiffusion in terms of the non-linear scaling 
in time of the mean squared displacement (MSD) 2] 



(r2(t)) ~ /i„t" with < a < 1, 



(1) 



where a is the anomalous diffusion exponent and Ka the 
generalized diffusivity of physical dimension cm^/sec". 
Subdiffusion ([1]) was reported for various microscopic 
tracers under the densely crowded conditions inside liv- 
ing cells [3- 7] and in control experimen ts M , [§| , as well 
as for proteins in cell membranes (13, Er These ex- 
periments demonstrate the ubiquitous presence of sub- 
diffusion of a large variety of tracers and crowded envi- 
ronments over many orders of magnitude in time, but see 
[l^ . Subdiffusion alters significantly the diffusion control 
of biochemical reactions, and its effects are therefore far- 
reaching for a wide range of molecular cellular processes 
[ist . While subdiffusion ([TJ slows down long-distance dif- 
fusional exchange and may affect surface-bulk exchange 
[3], it may indeed be beneficial for local interactions in 
cells d, [H^, [l^ . Depending on the magnitude of the ex- 
ponent a anomalous diffusion may effect the localization 
of objects such as chromosomes or membrane channels 
0, El, and impact on the formation and dynamics of 
membrane domains. 

Here we study in detail the diffusive behavior of lipids 
in bilayer systems through trajectory analysis from ex- 
tensive molecular dynamics simulations. We find that 
in all investigated bilayers the lipids exhibit subdiffu- 
sion up to a few nanoseconds, before a crossover to ei- 
ther normal diffusion or to persistent anomalous diffu- 
sion with a larger exponent. The observed behavior de- 
pends strongly on the phospholipid chemistry, their mix- 
ture with cholesterol, and the bilayer phase (liquid/gel). 
Subdiffusion ranges at least up to several hundreds of 
nanoseconds in the presence of cholesterols. Our analysis 



shows that the lipid motion is consistent with viscoelas- 
tic subdiffusion driven by correlated Gaussian noise in 
both liquid and gel phases, and thus provides a unified 
physical framework for lipid diffusion in membranes. 

Subdiffusion ([T]) is described by several prominent 
models based on different physical mechanisms 0, Il7| . 
fn continuous time random walks (CTRWs) jumps are 
separated by random waiting times r with heavy-tailed 
distributions ~ r-^"" [Hj. CTRW motion was 

identified for microbead motion in reconstituted actin 
networks [l^, lipid granules in cellular cytoplasm @, 
and protein channels in plasma membranes [f f]. In con- 
trast fractional Brownian motion (FBM) and fractional 
Langevin equation (FLE) produce ergodic subdiffusion 
^ with long-ranged anti-correlation (0 < a < 1) [1^ 



(Ar(t) • Ar(0)) - a{a - l)t"- 



(2) 



of spatial displacements Ar. FBM is defined by an over- 
damped Langevin equation driven by athermal, exter- 
nal Gaussian noise with power-law correlation. FLE is a 
generalized Langevin equation driven by the same noise. 
Due to its memory kernel the FLE describes thermal mo- 
tion of a particle in viscoelastic media [llj. While FLE 
in the overdamped limit produces FBM-like subdiffusive 
motion, below the momentum relaxation time, FLE mo- 
tion is ballistic. FBM/FLE motion describe subdiffusion 
in living cells of mRNA chromosomal loci [4|], and 
of lipid granules at longer times [ij, as well as the mo- 
tion of macromolecules in a crowded dextran solution Q . 
While sharing the scaling form of the MSD (H]), CTRW 
and FBM/FLE lead to completely different dynamics of 
diffusion-control. Knowledge of the subdiffusion mecha- 
nism in membranes is therefore vital to advance our un- 
derstanding of their physical and biochemical properties. 

Lipid bilayers are quasi two-dimensional, highly packed 
systems made up of phospholipid molecules, which un- 
dergo thermally driven lateral diffusion and thus con- 
stantly reorganize the membrane. The lateral MSD of 
membrane lipids typically spans three distinct regimes: 
short-time ballistic (a = 2), intermediate subdiffusive 
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FIG. 1: Final configurations of simulations of DOPC 338 K 
(Left), DSPC/cholesterols 338 K (Middle), and DSPC 310 
K (Right) in the liquid disordered, liquid ordered, and gel 
phases, respectively (note the difference in packing states). 
Each color represents a different phospholipid. Explicit wa- 
ter molecules correspond to the upper and lower transparent 
coatings. Cholesterols appear in white (Middle), Na"*" and 
Cl~ ions as blue spheres (Right). 
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FIG. 2: TA MSDs {S^{A)) of DSPC, SOPC, and DOPC in 
liquid phase b ilayers. Each curve represents the mean of indi- 
vidual 5^ (A) taken over all trajectories of 128 phospholipids 
in the bilayer. A. Cholesterol-free case. B. With cholesterol. 
The results were fitted separately to {S^{A)) = 4KoiA" in the 
regimes of short [0.01ns. . .Ins] and long times [10ns. . .100ns]. 
Fit results for a and Ka are summarized in Tab. |l] 



(0 < a < 1), and long-time Brownian motion {a = 1) 
[23', ^231 . The long-time diffusive motion of various kinds 
of phospholipid molecules in lipid bilayers has been exten- 
sively studied [25|, l26| . Diffusion of lipids in pure bilayers 
occurs both in the liquid disordered and the gel phases 
below the melting temperature, the latter with decreased 
diffusivity. Moreover, in bilayers mixed with cholesterols, 
the diffusivity of the lipids tends to decrease with higher 
cholesterol concentration. 

Lipid subdiffusion at shorter time scales is compara- 
tively poorly understood. In the traditional microscopic 
picture, the lateral movement of lipid molecules is as- 
sumed to occur through jumps when sufficient void space 
is thermally activated at nearest sites 2^ 27 1. Between 
jumps the molecule, caged by its neighbors, undergoes 
rattling motion. This CTRW-type jump-diffusion model 
has been used to estimate the diffusivities of lipids in the 
liquid-disordered phase and in bilayers containin g ch oles- 



terol [25|, |26|. However, atomistic simulations [28| and 
quasi-elastic neutron scattering experiment |29t showed 
that such jump-like displacements rarely occur, and the 
lipids move concertedly with their neighbors as loosely 
defined clusters. Moreover, conflicting results were re- 
ported on the stochastic nature of the lipid diffusion: 
Refs. [illsol demonstrated that the lipid motion is con- 
sistent with FLE dynamics, whereas Ref. [sif] claimed to 
observe CTRW-type motion governed by non-Gaussian 
fluctuations and scale-free rattling dynamics. 

Lipid bilayers of 128 phospholipid molecules were stud- 
ied by molecular dynamics simulations under periodic 
boundary conditions, for details see the Supplementary 
Material (SM) [s^. We used three pure single compo- 
nent lipid bilayers composed of DSPC, SOPC, and DOPC 
phospholipids in the liquid disordered phase [33[. We 



A pure membrane of 288 DSPC molecules was also simu- 
lated in the gel phase. Fig. [1] shows typical snapshots in 
the three phases. In this work we focus on the character- 
ization of the lipid diffusion. To that end we note that 
during the simulation the center of mass of the upper and 
lower lipid layers undergo free, independent translational 
motion (Fig. SI), as reported previously [3l|, |3J]. Free 



center of mass diffusion causes apparent normal diffusion 
of individual lipid molecules at longer times, irrespective 
of their actual diffusion characteristics. To avoid this we 
analyze the relative motion r{t) from the center of mass 
of lipids and cholesterols. Figs. S2, S5, and Sll in SM 
show sample trajectories. 

From individual trajectories r{t) we obtained the time- 
averaged (TA) MSD of lipids typically defined as 0,0, [13] 



52(A) 
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[r(i + A)-r(i)]2di, (3) 



also studied these systems with additional 32 cholesterols 
(20% molar concentration) in the liquid ordered phase. 



where A is the lag time and T the length of the trajec- 
tory (measurement time). Fig. [2] shows the mean {5'^ (A)) 
taken over the trajectories of all phospholipids, for the 
cases of DSPC, SOPC, and DOPC in absence and pres- 
ence of cholesterol. In each case, the result was fitted by 
((52(A)) = 4:KaA" at short and long times, respectively. 
The corresponding diffusion exponents a and diffusivi- 
ties Ka are summarized in Tab. HI In Fig. [21 the scaling 
behaviors for pure DSPC and DOPC at short (solid line) 
and long (dashed line) times are indicated. In absence of 
cholesterol, all three types of lipid molecules show simi- 
lar behavior: anomalous diffusion with exponent a ^ 0.6 
below a crossover time Tc ^ 10 ns, and normal Brownian 
motion beyond t^. The crossover time Tc roughly corre- 
sponds to the diffusion time of a lipid molecule needed to 
span its nearest-neighbor distance. The structural dif- 
ference in the tails of the lipids affects somewhat the 
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DSPC 


0.63(0.06%) 


0.032(2%) 


0.94(0.9%) 


0.020(5%) 


SOPC 


0.66(1%) 


0.038(2%) 


1.00(16%) 


0.020(60%) 


DOPC 


0.67(2%) 


0.043(2%) 


0.97(13%) 


0.028(16%) 


DSPC 


0.52(1%) 


0.013(2%) 


0.82(4%) 


0.0076(6%) 


SOPC 


0.58(0.4%) 


0.019(0.5%) 


0.87(5%) 


0.012(4%) 


DOPC 


0.61(0.5%) 


0.023(0.4%) 


1.04(22%) 


0.0098(41%) 



TABLE I: Exponent a and diffusivity Ka (nm'^/ns°) of lipids 
at short (s) and long (l) times. Statistical uncertainty in 
parenthesis. Boldface: systems containing 20% cholesterol. 



long-time diffusion, in particular, the values of . 

Fig.[33 shows that cholesterol significantly affects both 
the short and long time diffusion of the lipids. Especially 
for saturated DSPC with the smallest cross section area 
[ssj . we observe that below Tc w 10 ns a decreases to 
about 0.5 and a new subdiffusion regime emerges with 
a ^ 0.8 up to 100 ns. An additional 1 /is-long simulation 
confirms that this new regime above Tc is in fact a slow 
transition toward normal diffusion that lasts over hun- 
dreds of nanoseconds (Fig. SIO). In the displayed time 
window the behavior is well fitted by the scaling expo- 
nents ai listed in Tab. HI For unsaturated DOPC, the 
effect of cholesterols is small, albeit is significantly 
reduced, and no second subdiffusion regime occurs. 

The observed subdiffusive behavior of lipids is mainly 
attributed to the unique structural complexity of the lipid 
molecule. Spherical-shaped hard particle systems cannot 
have such a long subdiffusion regime and values of a as 
small as 0.5-0.6 (Fig. S18). To gain additional physical 
insight into the lipid motion we now check the detailed 
stochastic properties of the lipids. 

Time-averaged observables obtained from single trajec- 
tories provide information on the ergodic properties of a 
stochastic motion, and thus about the physical nature 
of the underlying dynamics. We call a process ergodic 
when the long time average of a quantity (e. g., t he MSD) 
equals the corresponding ensemble average Ullli]. For 
free CTRW subdiffusion [s^ ((52(A)) grows like A while 
the ensemble average (lit sc ales sub-linearly Free 
FLE motion is ergodicjii, , and (^2(A)j ~ A". The 
observation of a sublinear slope in Fig. [2] already indicates 
that the observed motion is not of CTRW type. This is 
further confirmed by the independence of 5"^ of the mea- 
surement time T, Fig.[3j3, in contrast to the T"~^ scaling 



of CTRW [17|, |38|. Fig. W shows the distrib ution (j) of 
trajectory-to-trajectory amplitude variations of ^^(A) for 
the 128 individual molecules of Fig.[3]4., as function of the 
dimensionless variable ^ — 5'^{/S)/ ((^^(A)). All curves are 
centered around the ergodic value ^ = 1. The broaden- 
ing of (j) with increasing A mirrors large fluctuations of 
5"^ at long A due to insufficient statistics when calculat- 
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FIG. 3: Stochastic a nalysis of DSPC in pure liquid phase bi- 
layer. A. TA MSD (52(A) of aU 128 DSPC molecules. B. 
(52(A) versus measurement time T, for A = 0.01, 0.1, 1 
ns. C. Normalized scatter distribution (^(^) of (52(A) versus 
C = (52(A)/((52(A)), for A=0.1 ns, 1 ns, 10 ns. D. Displace- 
ment autocorrelation function C stit) / C stiQ) of DSPC lipids, 
5t — 0.2 ns. The solid and dotted lines represent the fit- 
free forms of Cst{t)/CstiO) for FLE and CTRW (SM P]) at 
a = 0.63. E. Moment ratios {r* (t)) / {r^ (t))^ (regular) and 
{''max(i))/(''max(i))^ [mean maximal excursion (MME)] for 
DSPC, SOPC, and DOPC molecules. The horizontal line at 
1.49 distinguishes FLE motion ((r4ax(t)>/(»-nmx(i)>^ < 1-49) 
from CTRW ((rLx(t)>/(rLx(i)r > 1-49). The horizontal 
line at 2 is the expected value of the regular moment ratio for 
both FLE and CTRW motion. 



ing (5^. (j) also narrows at fixed A when T is increased 
(Fig. S3). These properties demonstrate that the lipid 
molecules perform ergodic motion different from CTRW. 
We obtained the displacement autocorrelation function 

Cstit) = ([r(i + 6t) - r(i)] • [r((5i) - riO)]) /St^ (4) 

for arbitrary time step 6t for several diffusion models in 
SM III [33. Normalized, Cstit) /Cst{0) is a fit free func- 
tion for given St and a. Fig. [3p shows Cst{t)/Cst{0) 
at St = 0.2 ns of DSPC lipids from simulations in the 
subdiffusion regime, alon g w ith theoretical results for 
CTRW and FLE motion [34|- We find excellent agree- 
ment with FLE motion. Here a — 0.63 was taken from 
the TA MSD. The lipid motion is thus anti-correlated, 
in line with Eq. ([2|). The behavior in Fig. [3jl) differs 
distinctly from free CTRW motion 17| where Cst{t) — 
for i > St (dashed line). Fig. [3jil also shows the 
moment ratios (r^ (t)) / {r^ [t))'^ and (r•;^ax(^))/('^max(^))^ 
where rniax(i) denotes the maximal distance of a given 
particle from its initial position reached up to time t 
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FIG. 4: DSPC molecul es in liquid phase bilayer mixed with 
cholesterols. A. (52(A) for 128 DSPC molecules (for choles- 
terols see Fig. S8). Average TA MSDs {S'^{A)) for cholesterols 
(blue circle) and DSPCs (red square) are shifted by a factor 
20. B. Distribution for of DSPCs and cholesterols. 



4l| . Moment ratios have unique values dependin g on 
the stochastic process |4l|, as summarized in SM |32j . 
Fig. ^ shows that {r"^ (t)) / {P (t))^ fluctuates around 
2, and (?'max(*))/(''max(^))^ decrcascs from « 2 to stay 
< 1.49, as predicted for FLE and violating CTRW. 

We conclude that the subdiffusive behavior shown 
above is robust, all analysis tools convincingly point to 
FLE motion. Analogous results were obtained for SOPC 
and DOPC molecules. The results are preserved at vary- 
ing temperature: temperature increase only leads to an 
increase of Kc^ in the Brownian regime (Fig. S4) . Consis- 
tent with previous studies [isl . |29| the collective motion 
of lipids exhibit a flow-like pattern (Fig. S16). 

As shown in Fig. [21 the diffusion of the lipid molecules 
is drastically changed by the presence of cholesterols. Is 
the stochastic nature also affected by cholesterols? We 
find that ergodicity of the motion is preserved, while 
cholesterols significantly affect the distribution 0(^) of 
5"^. Comparing with the pure bilayer (Fig. [3]C), with 
cholesterol 0(^) noticeably broadens (Fig. |4l3). Individ- 
ual lipids thus undergo considerable variations while dif- 
fusing in the presence of cholesterols, as seen in Fig. 
Moreover, with cholesterols the displacement autocorre- 
lation Cst{t) has a slightly deeper well, consistent with a 
stronger anti-correlation of the displacement (Fig. S6). 

For the motion of the cholesterol molecules themselves 
we find that their behavior is almost identical to the 
lipids', of FLE- type anomalous diffusion. Fig. |4j\ com- 
pares the TA MSD averaged over all trajectories {S^{A)) 
separately for DSPC and cholesterol molecules. While 
cholesterols universally diffuse faster than the lipids, their 
scaling behaviors are almost the same. Meanwhile, the 
scatter distributions of cholesterols are sharper than that 
of the lipids, implying that cholesterol diffusion is more 
uniform than that of lipids. This may be related to the 
fact that at any moment only some lipids are in direct 
contact with cholesterols which modifies their behavior 
from those of the remaining lipds. The displacement au- 
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FIG. 5: DSPC molecules in gel phase bilayer. A. {S^iA)). B. 
Displacement autocorrelation Cst{t)/Cst{0) with 6t = 0.2 ns. 
Solid line: theoretical result for FLE with a = 0.16. 



tocorrelation of cholesterols is hardly different from that 
of the lipids in the lipid-cholesterol bilayer (Fig. S9), and 
the moment ratios agree with FLE motion (Fig. S7). 

To obtain a full picture of the diffusive motion in 
lipid bilayers, we also studied the gel phase. For DSPC 
molecules we obtain: (i) Fig.[SJ^ shows that S'^ scales with 
as ~ 0.16 at short times and is thus remarkably smaller 
than the value as ~ 0.6 in the liquid phase. Moreover, 
in the gel phase the TA MSD remains subdiffusive with 
ai ss 0.59 beyond the crossover, (ii) The scatter distri- 
bution (j) of individual (5^ (A) shows that the lipid motion 
remains ergodic in the gel phase (Fig. S12) '42]. (iii) The 
results for the gel phase autocorrelation are consistent 
with FLE motion with exponent a « 0.16 (Fig. [5]B), as 
are the moment ratios (Fig. S13). (iv) Contrasting recent 
claims (sij . the rattling dynamics of lipids is consistent 
with FLE motion (Figs. S14, S15). 

In summary, we here report extensive molecular dy- 
namics simulations of lipid bilayer systems and the anal- 
ysis of individual trajectories using stochastic analysis 
tools. While we find a moderate dependence on the lipid 
chemistry, the effect of cholesterols is striking. Choles- 
terols effect more pronounced and persistent subdiffu- 
sion. FLE motion is identified as the unifying process 
for the motion of both phospholipids and cholesterols in 
liquid and gel phases. Our study thus provides an in- 
tegral picture of lateral motion of lipids by showing the 
compatibility of FLE-type stochastic motion of individ- 
ual molecules and their flow-like collective motion. 

Cholesterols significantly affect the phospholipids dif- 
fusion via increasing membrane packing and inducing 2D 
ordering [i^ (Fig. [1]). a is lowered significantly to ~ 0.5 
below Tc ~ 10 ns. In agreement, recent experimental 
and computational studies show that a decreases with 
increase of the concentration of proteins in the bilayer 
4J, |45| . Interestingly we observe a pronounced variation 
between individual lipid's motion, likely due to the asym- 
metric disturbance caused by cholesterols [i^ . While the 
slowing down of lipid diffusion by cholesterols is known 
from experiment j26|, [461] and simulations 3J] , the dra- 
matic effects of cholesterols on intermediate-time lipid 
diffusion have not been reported to our best knowledge. 

Given above results we speculate that in biomem- 
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branes, whose complexity is higher than the bilayers' 
studied here (e.g., larger number of lipid moieties, pro- 
teins, and higher cholesterol concentration), subdiflFusion 
may range to macroscopic times, thus altering our cur- 
rent view of membrane dynamics. Single particle track- 
ing together with advanced simulations techniques and 
stochastic analysis tools are promising methods to ex- 
plore this intriguing possibility. 
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